Abstract 



Kinetics of the first order phase transition is investigated. In some 
rather wide spread situations the property of avalanche growth of the ob- 
jects of a new phase doesn't take place. This radically complicates the 
theoretical description of the nucleation kinetics and this case was com- 
pletely investigated here. The power-like dependencies with decreasing 
absolute value of derivative on time were chosen as the model laws of em- 
bryos growth. All main characteristics of the nucleation period including 
the spectrum of sizes of the new phase objects were obtained. An analyt- 
ical approximate description of nucleation in different external conditions 
was fulfilled. The errors of such description are found and proved to be 
small. 
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1 Introduction 

Kinetic description of a system under the phase transition is typical namely for 
the phase transition of the first order. The driving force of the first order phase 
transition implies the existence of metastability of a mother phase. Namely 
this metastability can be regarded as a "source" of the substance transforma- 
tion from an unstable "mother" phase to a stable "new" phase. Thus, kinetic 
investigation of the phase transformation is rather actual for the phase tran- 
sition of the first order. For example, such an important characteristic as the 
number of objects of a new phase can be calculated only on the base of kinetic 
consideration. 

Traditionally the physical model situation for the first order phase transition 
was the case of the supersaturated vapor condensation into the phase of liquid 
droplets. It was chosen as the base for theoretical and experimental investigation 
due to the rather simple description of the bulk phase. Namely in this case 
the classical theory of nucleation was created in [Q. This gave the base for 
construction of the kinetic description of the global evolution of the system 
under the first order phase transition which was done in 0], That's why we 
shall borrow the terminology from the case of condensation @], 0, For example, 
the supercritical embryos of a new phase will be called as droplets, the process of 
the phase transformation will be called as condensation, the period of formation 
of the main quantity of the droplets will be called as nucleation, etc. 

One of the main features which allows to fulfill analytical constructions in 
0], |S] was the property of the avalanche consumption of the mother phase by 
droplets. This means that the quantity of the mother phase molecules absorbed 
by a small droplet is many times smaller than the number of molecules absorbed 
during the same time by a droplet of a big size. 

In a free molecular regime of the substance exchange one can write that 
absorption coefficient W + is proportional to the surface of droplet, i.e. to the 
square of the droplet radius R. Then W + grows in time t under the constant 
power of metastability as W + ~ t 2 and the number of molecules v accumulated 
by a droplet grows as t 3 . The big power 3 allows to apply in g|, p] some effective 
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methods which can be reduced to the first steps in some iteration procedures, 
formulated there. 

Under the diffusion regime one can get v ~ i 3 / 2 and even this not so rapid 
growth is also sufficient for the application of the methods from 0] based on the 
first iterations. 

All possible regimes of the substance consumption lies between the free 
molecular regime and the diffusion regime and namely this property allows to 
establish in |S] the general feature - the avalanche consumption of metastable 
phase. This property isn't only the key to get the numerical results but allows 
to combine the collective regime of the metastable phase consumption and the 
existence of the density profiles in a unique approach [3], [S]. This possibility 
plays the fundamental role in the justification of approaches formulated in 

But sometimes the rate of droplets growth decreases in time. This property 
doesn't allow to apply iteration methods from 0], 0- This situation requires a 
special consideration which will be given below. The situation with decreasing 
rate of growth takes place rather often in the processes of cementation, structural 
transformations, morphological transitions, etc. 

Rigorously speaking one has to state that the avalanche consumption is 
the characteristic feature of phase transition, the justification is based on the 
geometrical features (i.e. on the dimensionality of space). So, the situation 
with the decreasing rate of growth isn't characteristic for the pure transition in 
a rigorous sense. But in practice the decreasing rate of growth is also considered 
in frames of the first order phase transition in an ordinary sense. So, it is worth 
doing here through a presentation of the analytical theory of the process. 

We shall consider two types of external conditions which seems to be typical 
and are widely spread both in theoretical and practical investigations. The first 
type is a situation of decay. It takes place when the power of metastability is 
attained in the system rather rapidly and later there is absolutely no external 
influence upon the system. The power of metastability is fully determined here 
by the condensation processes. 

The second type typical external conditions are conditions of dynamic type. 
Really, ordinary practically all external conditions have rather smooth behaviour 
in time. It means that the metastability is created gradually. Formation of 
droplets and mainly the process of growth consume vapor and the power of 
metastability begins to fall. As the result the process nucleation stops and one 
have the formed spectrum of sizes of droplets with rather simple law of evolution. 

Now we shall discuss the law of the droplets growth. Here we shall use the 
power-like dependencies of the droplets growth. The power-like dependencies 
are explained by the absence of some characteristic dimensional scale like in the 
scaling theory for the second order phase transitions. The diffusion regime and 
the free molecular regime have also the power-like laws of growth. So, we shall 
also take the power-like regimes here. Namely 

v = t a 

with parameter a. We shall call the power dependencies with < a < 1 as the 
slow rates of growth. 
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We shall use all results and definitions formulated in [3J, 0], M- But, cer- 
tainly the methods used in these papers have to be reconsidered. They can not 
be directly applied under the slow rate of growth and a new methods have to 
be developed. 

We investigate the case of homogeneous phase transition. The heteroge- 
neous case can be investigated analogously to 0]. In 0] it is shown that the 
homogeneous case can be effectively applied as the base for the consideration of 
heterogeneous case. The same takes place here also. 

2 Evolution equations 

The situation with the slow rate of growth is more complex because now not only 
the droplets with relatively big sizes are the main consumers of vapor. Here all 
droplets takes place in the consumption of vapor. But there is one simplifying 
feature which states that now there is no need to consider the density profiles or 
profiles of metastability [3J, 0. Really the characteristic scale of the diffusion 
blurring is 

Rdiff ~ vt — V 

where t is the current moment of time, t' is the time of the perturbation ap- 
pearance. The same is valid for the thermal blurring. So, when a < 3/2 the 
diffusion blurring destroys the density profile and we come to the collective 
regime of vapor consumption (see [S])- 

The procedure of derivation of the evolution equation is absolutely the same 
as in 0], one has only to substitute the power 3 or 3/2 by the new power a. The 
statement about the leading role of the supercritical embryos in the substance 
balance used in 0] is so strong, that it remains valid until a — > 0. But the case 
a —* can be investigated explicitly without this assumption which will be done 
separately (due to the simple conservation laws) below. The same is valid for 
the statement about the applicability of the quasistationary approximation for 
the rate of nucleation (the rate of appearance of droplets). 

The power of metastability will be characterized by the value of supersatu- 
ration 

C = (n- n^/rioo 

where n is the molecular number density of the mother phase, is the same 
value for the saturated mother phase. 

The evolution equation for the situation of decay can be written as 

G = A f (z-x) a exp(-G)dx 
J o 

with a positive parameter A, which has the sense of the spectrum amplitude. 
The spectrum looks like 

/ ~ exp(-G) 
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In the situation of dynamic conditions we can linearize the power of metasta- 
bility in the absence of nucleation as we have done in £0 . Then we have 

G = A { {z-x) a exp(bx - G)dx 

J — oc 

with the additional parameter of the mentioned linearization b. Here 

/ ~ exp(bx — G) 

In the derivation of these equations we have used the following approximation 
for the stationary rate of nucleation 

/ s (0=/,(Co)ex P (£>(C-Co)) (1) 

as a function of the supersaturation £ where index characterizes the base of 
decomposition. Here D is the derivative of the free energy of a critical embryo 
formation over £. 

This approximation isn't something special but simply the Klapeiron-Klausius 
equation. Ordinary this equation is used for the dependence of over tem- 
perature but there is no difference to what characteristic we shall apply this 
procedure of approximate derivation. We have to mention that renormaliza- 
tions to describe the nucleation kinetics in non isothermal conditions are based 
on the Klapeiron-Klausius equation for the density of the saturated vapor. Here 
the analogous renormalizations in order to take into account the non isothermal 
effects can be also done. 



3 Decay of metastable state 

At first we shall investigate the evolution of the system under external conditions 
of the decay type. Certainly, one can act in the manner analogous to the 
approach of the universal solution [Sj for the case of the avalanche consumption 
of the mother phase. After the evident rescaling z — > A 1 / i z, x — > A 1 / i x one can 
come to 

G= [ (z - x) a cxp{-G)dx 
Jo 

with no parameters. Then G is the universal function and the zero momentum 
Mo = Jq is the universal constant. Since M is the coefficient in the asymptot- 
ically leading term Mqz" in expression for G after the end of nucleation we have 
got all information necessary for the further evolution. This solves the problem 
of the nucleation description. 

Having determined the duration of the nucleation period Sz as 

N(Sz) = AT(oc)(l - e) 

with some small e and N(z) as the number of droplets appeared before z one 
can get for Sz some fixed constant. The value N(z) ~ / exp(— G)dx is also the 
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universal function and Sz as the root of equation with no parameters will be the 
universal function. 

Now we shall determine the form of the size spectrum of droplets explic- 
itly. The iteration procedure isn't effective because the initial approximation is 
too rough. To propose another iteration procedure we shall establish the new 
property of the droplets growth typical for the regime of th e slow growth. 

The property of effective size of growth 

Let us consider the law of the the vapor consumption by a separate droplet. 
It is shown in fig. 1. For every moment z (or t(z)) we can find the boundary zq 
with two properties 

• zo <C z 

• v{z) — v(z ) <C v{z) 

One can see that for z — az /p, P ~ 2 -r- 3 both properties are satisfied. 

It means that the qualitative picture is the following: rather soon the droplet 
accumulates the main quantity of vapor and during the rest of the evolution time 
the accumulation isn't too important. Certainly, this moment of accumulation 
grows in time. 

This picture allows to suggest the following initial approximation: 

For every z one can imagine that practically immediately the number of 
molecules of a given droplet (with coordinate z) attains the value v(zq) and then 
the droplet doesn't grow. 

One has to stress that for different z the values zo are also different. 

This property allows to use the approximation with a = as initial approx- 
imation, which has an evident advantage. 

Here one can see the property of a real collective consumption of vapor - 

Practically all droplets ( except small ones with sizes less than z$ ) consume 
vapor in approximate equal quantities. This is really the collective consumption 
of vapor (Earlier this term was used simply for vapor consumption in a fixed 
point by many droplets in different quantities. Now droplets consume vapor in 
equal quantities and we shall speak about the equal collective consumption.) 

The case a = can be solved analytically, because it can be reduced to the 
ordinary differential equation of the first order. This leads to 

exp(- G ) = ^ 

One has to note that the last equation leads to the infinite number of droplets 
appeared in the process of nucleation. This is an error. This error is initiated 
by the inapplicability of approximation Jy| at small supersaturations. Really, 
at G ~ n — tt-oo the rate of nucleation has to be vanished, but approximation 
Q gives the finite value of the nucleation rate at any G. 

Nevertheless we aren't interested at the long tail ~ z^ 1 of the size spectrum 
because this tail can not be directly seen in experiment and all integral values 
can be explicitly calculated. Really, the total number of droplets can be easily 
obtained, it is convenient to do even in initial variables. 
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The case a = corresponds to a fixed number Vfi n of a number of molecules 
inside the droplet. This quantity doesn't depend on time, in rescaled units 
Vfin = 1- The total number of droplets is 



,, _n-n OQ 

N total — 

Vfin 

Formally we can choose the cut-off of the spectrum at the size where N{z) 
equals to N to tai- 

Iteration procedure 

One can see that even in the case of small positive a the spectrum of sizes 
is well localized which allows to use the iteration procedure 



G i+1 = (z- x) a exp(-Gi)dx 
Jo 



Here the lower index denotes the number of iteration. As initial approximation 
we choose here the analytical solution for a = , i.e. 

G = - ln(x + 1) 

This radically differs from consideration of the situations with avalanche regimes 
of growth, where Go = 0. 

The first iteration can be easily calculated. We have 

f z 1 
Gi = / (z-x) a dx 



(i 



1 + x 



which can be easily calculated 



z Z 

G t = -__$(___, l,a+l) 

2+1 Z+l 

where 

oo 



z k 



is the standard special function 

The spectrum based on G\ , i.e. exp(— Gi) can be considered as a rather 
precise approximation. It can be shown analytically and numerically. Namely, 
the relative error in the number of droplets, i.e. in 



N(z)= / exp(-G{x))dx 
Jo 

is less than 0.05 for a < 0.5. The values a from 1/2 up to 3/2 can be con- 
sidered analogously on the base of iterations started from solution with a = 1. 
This solution can be also found analytically since the evolution equation can be 
reduced to the second order ordinary differential equation with no explicit de- 
pendence on the argument, which can be easily integrated. Then is is necessary 
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only to take one iteration step and then to integrate exp(— G\) over time to get 
the number of droplets. 

Figure 2 shows the relative error in the total number of droplets as a function 
of a. The integration is fulfilled up to z = 10 which is absolutely sufficient for 
any case except a — and in this case the error is zero (certainly with artificial 
cut-off the same will be numerically). The deviation from zero at small a is 
initiated only by a finite step of integration. One can see that the relative error 
is rather small. This can be shown not only numerically but also analytically. 

Figure 3 shows the characteristic forms of spectrums. Two cases are shown - 
a = 0.2 and a = 1. One can see two pairs of curves coinciding at the beginning 
and approximately at the end. Every pair corresponds to one case and presents 
the precise solution and our approximation. The case a = 1 has analytical 
solution and it is out of consideration here but one can see that even in this 
case our approximation works good. In this case the back side of spectrum is 
more sharp than for a = 0.2. One can also notice the relative invariancy of the 
form of spectrum to the concrete choice of a for small a. 

One has to mention that the consumption of vapor by all droplets requires to 
reconsider slightly the characteristic size of consumers in the statements about 
the quasistationarity of the nucleation rate and about the size of the main 
consumers of vapor [3] . But even with these modifications both statements can 
be proven (and under dynamic conditions also) 

4 Dynamic conditions 

The simple rescaling z — ► zA 1 / 01 , x — ► xA x l a brings evolution equation to 

G= f (z - x) a exp(bx - G)dx 

and allows to cancel parameter A. But it isn't too simple to cancel parameter 
b because earlier [5] the condition, which allows to cancel this parameter was 
the following: The point of decompositions has to be chosen as the maximum 
of supersaturation (or as the maximum of /) 

In situation with a = one can not satisfy this condition. Moreover, / 
increases in time. The analytical solution in the case a = is the following 

/ = exp(6x-ln(^PM + l))^/ (2) 

This solution can be got after differentiation of the evolution equation which 
brings it to 

— — = cxp(6a; — G{x)) 
ax 

One can easily integrate this equation 



exp(G)dG = J cxp(bx)dx 
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and come to (J2J. 

When a =/= one can see the maximum of / and when a = the explicit 
analytical solution has been presented. 

For a > one can put condition on maximum. Then for a^O we have the 
universal solution depended on a. This solution satisfy the following universal 
equation 



o 

a-l 



-ln(/) + a / (z - x) a ~ 1 f{x)dxz = I {z-x) a f(x)dx 



the ideology here is absolutely analogous to the already investigated situations. 
Iteration steps 

Now we shall present the methods to get the explicit form of the size spec- 
trum. The main feature which has to be taken into account is the long infinite 
tail of the spectrum for a = 0. The explicit form of spectrum in this case for 
6 ~ 1 is drawn in figure 4. One can see that the curve is slightly blurred because 
here an approximation which will be discussed later is also drawn. 

The most interesting region here is the region of the rapid increasing of /. 
If we choose b = 1 we get this region at x ~ 0. Later we shall see that it is 
reasonable to put b ~ 1 for all a. 

One can formulate an important feature: Certainly, for a > the spectrum 
/ lies lower than the spectrum /o in the case a = 0. Having used the property 
of effective size of growth one can come to conclusion that the maximum of / 
lies near x = when 6=1. So, then one can put 6=1 for all cases with small 
a. 

We shall use /o as initial approximation in the iteration procedure 
G i+ \ = (z — x) a exp(6x — Gi)dx 

J —oo 

Go = - ln(/ ) + bx 
The first iteration can be presented as 

d= T (z-x) Q (^M + l)- 1 exp(6x)dx 

J — oo ^ 

To calculate G\ one can act in two manners. The first one is to invent approxi- 
mation for (z — x) a . as we have already mentioned this function has behaviour 
drawn in figure 1. We see that that it can be well approximated by a straight line 
as it is drawn in figure 5. This approximation allows the analytical calculation 
of G\ since the integral can be reduced to 

-dx 



and to 



(a; + ci) 2 - c 2 

ln(x + c 3 ) 1 c 5 x 2 2 

-— — -2 dx == — arctan In c 3 - c 4 + c 5 

(X + C4) z - C5 2C5 C| + ci, + C4X 
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-^-Cl 2 (2e + 2$) + ^-Cl 2 (2e + 2$) - -^-Ch{ir - 29) + ^Ch{ir - 29o) 

ZC5 ^£-5 

Here 

tan 9 = 

tan9o = — 

C5 

, c 3 - c 4 
tan $ = 

C5 

and 



C7 2 (z) = - f ln(2sin(^)da 



for x < and 



is the Klauscn integral which is the standard special function. 

The second method which suggests an approximation for /o seems to be 
more attractive. This approximation is the following 

/o = z ; 7 r ~ fap 
1 + exp(— x) 

where ^ 

f ap = exp(x) - exp(2x) + - cxp(3x) 

fap = 1 - exp(-a;) + - cxp(-2a;) 

for x > 0. 

Namely this approximation together with precise solution is drawn in figure 
4. Certainly at x — there is no difference between precise result and ap- 
proximate one. With the help of this approximation all integral terms can be 
calculated in terms of Gamma-function. 

When the first approximation G\ is calculated we see that the specrtum f\ = 
exp(a; — G\) in the first approximation has one maximum f\ max at x± m ax and 
at the region x < x\ max the spectrum /1 approximates the real spectrum quite 
satisfactory. Moreover we have x max w x\ max, where x max is the coordinate 
of real maximum and f m ax ~ fi max where fmax is the amplitude of a real 
spectrum. These facts can be proven analytically. The form of size spectrums 
for characteristic values of a are drawn in figure 6. 

Gne can see one important feature here. The spectrum in the first iteration 
rapidly turns to zero while the real spectrum has very long and smooth tail. This 
occurs because the amplitude of the spectrum (and also of the first iteration) is 
seriously less than the amplitude of the solution for a = 0. Now we can correct 
this error. When a is growing then x 1 ma x goes away from x max ; f± max goes 
away form f max and fo(xi max ) goes away from f\ max . But even when a = 0.5 
(this is the realistic boundary of our constructions because for bigger a one can 
use decompositions starting from a = 1) the difference between fo(xi max) and 
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fi max is small. This allows to suggest the cut-off of the zero approximation 
fo cut which is 

fo cut = fo 

for x < xi max and 

fo cut — fo(%l max) fl max 

for X > X\ max- 

On the base of this initial approximation one can reproduce all constructions 
and get the spectrum fx in the first approximation, which lies higher than the 
"previous first approximation". It is drawn in figure 7 where this curve lies 
higher then the previous first iteration. 

One can analytically show that the first advanced iteration describes the 
form of the hill (the head of the spectrum) satisfactory. But there remains the 
long tail and it is absolutely impossible to describe t on the base of the iteration 
method because the duration of the back tail can be infinitely long. Then one 
can propose come asymptotic methods to solve this problem. 

Asymptotics 

We shall construct solution for x » 1. At first we can see that it is possible 
to write approximately evolution equation as 



G{z)= / (z-x) a exp(x-G(x))dx 

instead of 



G(z) = [ (z- x) a exp(x - G{x))dx 

J — CO 



Now we shall use the property of effective size of growth. Since the amplitude 
of spectrum doesn't vary too rapid during this time one can take away the size 
(z — x) a out of the integral and speak about the mean size p. Then 

G(z) —pi exp(x — G(x))dx 
Jo 

For p one gets the following expression 

__ Jo(z-xrdx _ z» _ a 
J Z dx a + 1 

Having differentiated G over z one has to notice that p — z a is rather slow 
function and there is no need to differentiate it. Then 

^ = z a ^- f cxp(x — G(x))dx = z a cxp(z — G(z)) 
dt dz Jo 

After intergation we have 

j cxp(G)dG — p J cxp(x)dx + const 
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and 

G = ln(l + pexp(z)) 
This solution has asymptotics 

G -> Mp) + z 

This asymptotics leads to the following expression for the size spectrum 
/ <~ exp(a; — G) = exp(z — ln(p) — z) = 

and finally 

Here one can introduce the arbitrary shift 7 and the arbitrary amplitude (3. 
Then 

/ V— 

(z - 7)° 

These parameters can be determined by requirement of the smoothness of spec- 
trum and its derivative in the point of transition from the advanced first iteration 
to the asymptotics. Were shall we choose this point? 

At z = the spectrum begins to grow, at z = z max it attains the maximum, 
it is reasonable to imagine the hill as symmetric one an dto say that at z = 2z max 
the hill is over. Namely at this point we can speak about the beginning of the 
asymptotic. 

The conditions 

A P 



dz (z - 7)"+! 

and 

f P 
h 



at z = 2z max give the following expressions 

7 = 2z max + afi(z max )/(dfi/dz)\ z = Zmax 

(3 = h{Zmax){z ~ lT 

The characteristic forms of spectrums and analytical approximations are 
drawn in figure 8. The letters A-E denote the pairs of curves corresponding to 
different a from 0.1 to 0.5 with a step 0.1. The capital letters mark the precise 
numerical solutions, they lies near approximate solutions. 

The accuracy of the theory can be estimated by the error in the droplets 
number 

\N-N ap \ 
6 N 

is drawn. Here N is precise value found from numerical solution and N ap is the 
number of droplets found from the presented approximate solution. 
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In figure 9 the maximum of e over z is drawn as a function of a. It can be 
seen that it is rather small. Here the maximum value of z is chosen as 30. 

Now we shall analyze the accuracy of calculations. In figure 10 the relative 
error for N(z) at some final value is drawn. Here we consider the "final" values 
for N. Certainly, the terms "the final values" are illegal, because it is clear that 
for a < 1 the finite number of droplets can not satisfy the balance of substance. 
So, we need to throw away th e tail, which is thin but still infinite (it contains 
infinite number of droplets). 

In consideration of extremely long asymptotic tails we have to take care 
about the accurate behavior of asymptotics. The combination of iteration solu- 
tion and asymptotics has to be performed here in a slightly other special point. 
We shall mark this point as Zb OU nd- 

To get Zbound we shall fix the beginning of nucleation more precise. On the 
base of iterations we get fmax- Then we get the value f s t of the amplitude 
at the beginning of nucleation as f s t = fmax / exp(l). Then we can get the 
coordinate z s t of the beginning of nucleation as z st = ln(/ st ). Then Zbound will 
be calculated according to the old recipe but with a new time of the beginning 
of nucleation 



In figure 11 the values of the droplets number are compared at a "final" 
value Zfi n . Two values z/,„ = 50 and Zfi n — 70 are chosen here. The value 
of e is drawn. One can see three broken lines. The shapes of two of them are 
approximately the same at small a - these lines are the relative errors in the 
droplets number for different Zf in . It is clear that the dependence on Zf in at 
such values is very smooth and we practically attains the limit case (precisely 
speaking, it can not be done). The step in calculations dx was chosen as dx = 
0.05. The third line at small a lies below and has a more smooth shape. It is 
the same value e for Zfi n = 50 but calculated with a step dx = 0.025. It isn't 
too far from the two mentioned curves. So, the necessary accuracy is attained. 

Now one can analyze the behavior of the error. The decrease of the error 
ata« 0.7 is caused simply by compensation of different sources of errors, the 
value \N — N ap \ changes a sign here. More interesting fact is the reduction of the 
error at a ~ 1. The reasons are the following: Every spectrum has the tail and 
the head. The head is rather short and it is localized well and can be described 
with the help of iteration method. The tail has to described by asymptotics. 
The main source of error is th absence of account of the influence of surplus 
substance appeared in the head of spectrum in evolution at asymptotics. When 
we are going to extremely long tails this influence will disappear. 

When a attains 1 the quantity of the droplets in tail at intermediate z is small 
in comparison with the quantity of droplets in the head. So, it is reasonable 
to take in to account only the head with the help of iterations (Certainly, now 
the number of droplets wil be counted on the base of iterations not only until 
Zbound but also fo iz > z bound ). So, now 



Abound z st ~t~ 2(Z max 



z — st 
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fl = cxp(.T - gi) 
9i = I x) a fo{x) 

J — OO 

f = max fx 



for x > x m ax i and 



fo = f 



for x < Xmax 1 ■ Here maxfi is the maximal value of f\ and x m ax 1 is the 
corresponding argument. To get these values it is necessary to calculate 

fl = exp(a; - g x ) 
.91== f {z-xYMx) 



where /o is the analytical solution corresponding to a = 0. 

The result is drawn in figure 12. Here we draw the relative error e at z = 50. 
The step of calculations was dz — 0.05. We see that even for a = 3/2 the 
result is good even without asymptotics and even with initial approximation 
corresponding to a = 0. For a > 3/2 the result was presented in [o]. All 
situations are solved now. 

We have to stress that initial approximation as the analytical solution for the 
case a = 1 is better than initial approximation corresponding to a — 0. Then 
the accuracy will be also better. Now we shall present the analytical solution 
for a = 1. The evolution equation is 



(z — x) exp(x — g(x))dx 
For <f> = —x + g this equation is 

(z — x) exp(—4>(x))dx 



Having differentiated two times we get 

d 2 6 



-—p = exp(,-0j 

This second order differential equation doesn't contain the argument explicitly 
which allows to integrate it. Let <fi be the argument u, d<p/dz be the function y. 
Then d 2 4>/dz 2 — ydy/du. Integration gives 

y 2 

-— = exp(u) + ci 
From the boundary conditions we get 

1 

Cl = ~2 
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Then 

d(j)/dz = \fc\ — 2cxp(0) 

Integration gives 

d<p 



/ 



sjcx - 2exp(» 



This integral can be easily taken analytically which gives the analytical expres- 
sion for the spectrum in this case. The constant C2 can be got from the limit 
behavior (f> — > — x when x — > — oo. 

This solution is very fruitful for description of situations with a close to 1. 
Results are given in figure 13 which is analogous to the figure 12 but with another 
initial approximation. One can see here the curve A, which demonstrates the 
error of the first iteration in direct iteration method. Here 



/oo 
fi(x)dx 
-oo 

fi = exp(a;-5i) 
.91 == f (z - xTfo(x) 

J — oo 



and /o is the analytical solution corresponding to a = 1 . 

We see that the error is even greater than the error in figure 12. It is so 
because simply the size (the amplitude, not the shape) seriously diminishes. 
This is the same as we have seen for smaller powers. This diminishing is the 
main source of error. So, we need to reexamine the initial approximation. 

The first way is to act in a style like it has been done with the advanced 
iterations quite above. Here 



/oo 
fi(x)dx 
-OO 

fi = exp(x-gi) 
91 =/'(*- xTfo(x) 

J —oo 

fo = maxfi 



for x > x max i and 



fo - h 

for x < x max i. Here maxfi is the maximal value of /i and x max i is the 
corresponding argument. Then 

/i = exp(a; - g x ) 
9i = I (z~ x) a f (x) 

J — OO 
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and fo is the analytical solution corresponding to a = 1. The results are given 
by the curve B. 

We can act also in another manner. We can rescale the number of molecules 
in a liquid phase g to have the approximate equal amplitudes. One can calculate 
constants 



and instead of g consider g/q (both for the current a and for a = 1. This 
corresponds to the approximate equal g at z = 0. The error is drawn as the 
curve C. One can see that the error seriously diminished. 

One can also require the approximate equality in derivatives of g on x at 
t* as it is considered in the balance for establishing of i* (see @]). Then the 
constants q will be 



The result is drawn by the curve D. It is seen that the error is practically the 
same. It corresponds to the property of approximate universality observed in 
PJ. Now we see that this approximate universality goes also for the case of 
intermediate 1 < a < 3/2. 

Here and above the values of q were calculated on the base of the ideal 
supersaturation, i.e. on the base of exp(x). The results will be even better if 
we take q calculated of the base of solution at a = 1. 

One can not calculate here iterations analytically, but can act in a manner 
presented in [3]. Certainly, the leading term here will be zN where z is the 
coordinate of the maximum of peak of spectrum. 

The situation of decay is much more simple and we can use the standard 
iteration solution given by formulas 



where fo is the analytical solution corresponding to a — 1. The error is drawn 
in figure 14. It is small. The slight diminution when a grows is caused only by 
numerical errors. 

One has also to keep in mined that in situation with long tails one has to 
take into account the higher derivatives of the free energy of the critical embryo. 
This account is rather easy since we know the general solution. It can be done 
both by iteration approach and by perturbation technique. 

Also one has to take care about the validity of linearization of ideal super- 
saturation for the long tails. But it very simple to correct the asymptotics for 
deviations from the linear case. 






h 



exp(-si) 
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Figure 1 

Linear approximation of the size of growth 
Power a = 1/8 

The bold line presents the approximation 
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Figure 2 
Relative error 
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/ = exp(-G) 
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Figure 3 

Forms of spectrums 
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Figure 4 

Forms of spectrum and approximation 
for a = 
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Figure 6 

Forms of spectrums in the first iteration 
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Figure 7 

Forms of spectrums in the first modified iteration 
a = 0.2 
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Figure 8 

Forms of spectrums 
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Figure 3 

Maximal relative errors 
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Figure 

Relative errors 
The lenght is 70 
The step is 0.05 
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Figure 

Relative errors for pure iterations 
The lcnght is 50 
The step is 0.05 
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Figure 

Relative errors for pure iterations 
The lenght is 50 
The step is 0.05 
Basic power 1 
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Figure 

Relative errors for pure iterations 

The lcnght is 50 

The step is 0.05 

Basic power 1 

decay 
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